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O 



Motivated by the work of Guillot (2010), we present a semi-analytical formalism for calcu- 
lating the temperature-pressure profiles in hot Jovian atmospheres which includes the effects 
of clouds/hazes and collision- induced absorption. Using the dual-band approximation, we as- 
CO ' sume that stellar irradiation and thermal emission from the hot Jupiter occur at distinct wave- 

lengths ("shortwave" versus "longwave"). For a purely absorbing cloud/haze, we demonstrate 
Q, ' its dual effect of cooling and warming the upper and lower atmosphere, respectively, which 

rj"] I modifies, in a non-trivial manner, the condition for whether a temperature inversion is present 

^ ' in the upper atmosphere. The warming effect becomes more pronounced as the cloud/haze 

^ I deck resides at greater depths. If it sits below the shortwave photosphere, the warming ef- 

fect becomes either more subdued or ceases altogether. If shortwave scattering is present, its 
dual effect is to warm and cool the upper and lower atmosphere, respectively, thus counter- 
acting the effects of enhanced longwave absorption by the cloud/haze. We make a tentative 
j^ I comparison of a 4-parameter model to the temperature-pressure data points inferred from the 

observations of HD 189733b and estimate that its Bond albedo is approximately 10%. Besides 
their utility in developing physical intuition, our semi-analytical models are a guide for the 
^ i parameter space exploration of hot Jovian atmospheres via three-dimensional simulations of 

d^ ' atmospheric circulation. 

OS '. 

^^ , Key words: planets and satellites; atmospheres - radiative transfer 
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1 INTRODUCTION 



The study of extrasolar planets has swiftly tran sition ed from discovery to characterization. Following the pioneering detection of the 
atmosphere of a hot Jupiter by Charbonneau et ajj 1 20021). resea rchers have used the tran sit method to infer the presence of clouds and 



hazes in these atmospheres fe.... lPont et alJ2008l : l7inT:ralJ201lL rSee lportnev et alE(^ and references there.n for a detailed comparison 



of the observed atmospheric properties of hot Jupiters to spectral models.) Such discoveries motivate the improvement of atmospheric 
j^ models to include the effects of clouds and hazes. In particular, the simplicity of one- or two-dimensional models allows one to develop 
physical intuition prior to including new physical effects in (expensive) t hree-dimensional simulations o f atmo s pheric circulatio n. The main 
purpo se of the present study is to generalize the formalism presented in lHubenv, Burrows & Sudarskvl j2003h . lHansenl yOOSJ) and lGuilloli 
( 120100 to include several new effects: scattering, collision-induced absorption and additional sources of longwave absorption. The joint 
consideration of some of these improvements allows us to understand the effects of clouds/hazes on the temperature-pressure profiles of 
hot Jovian atmospheres. The simplicity and versatility of our semi-analytical models allow for an easy comparison to present and future 
observations aiming to characterize the atmospheres of hot Jupiters. 

Absorption and scattering introduce competing effects. Enhanced longwave absorption by the cloud/haze has the dual effect of cooling 
and warming the upper and lower atmosphere, respectively. By contrast, if the cloud/haze scatters in the shortwave, it has the counteracting, 
dual effect of respectively warming and cooling the upper and lower atmosphere. All of these effects combine to make the generalization of 
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Figure 1. Dependence of the Bond albedo A. and the coefficient in Tii-i- (or Tcq), (1 — ^)^''*, on the scattering parameter ^. 

the condition for the presence or absence of a temperature inversion difficult, as opposed to its crisp one-parameter description in a cloud-free 
scenario (Hubenv, Burrows & Sudarskv 2003; Hansen 2008; Guillot 2010). As a demonstration of the utihty of our models, we tentatively 
compare our semi-analytical, global-mean temperature-pressure profiles to the data points inferred from the observations of HD 189733b 
and estimate that its Bond albedo is about 0.1 (Figure[8j. 

In ij2] we describe the derivation of a general equation for the temperature-pressure profile in a hot Jupiter atmosphere, which allows for 
an arbitrary functional form for the longwave opacity kl . In ij3] we derive more specific forms of the equation assuming various functional 
forms for kl . Examples of temperature-pressure profiles are calculated in iJ4] We discuss the implications of our results in ij5] Table[T]contains 
a list of the commonly used symbols in this study. Our main technical results are stated in equations l l45t and l l50t . Appendices IaI IbI Icl and 
|D]archive technical information relevant to the formalism presented in !j2l 



2 GENERAL FORMALISM FOR TEMPERATURE-PRESSURE PROFILE 

In this section, we aim to analytically derive a general expression for the global- mean temperature-pressure profile of a hot Jupiter 
atmosphere by generalizing the formalism presented in Section 3.2 of iGuilloJ l[201O). We assume that stellar irradiation ("shortwave"; 
denoted by "S") and thermal emission ("longwave"; denoted by "L") from the hot Jupiter peak at distinct wavelengths, i.e., the dual-band 
approximation. The zeroth, first and second moments of the specific intensity are represented by J,^, H^ and K^, respectively, where v 
denotes the photon frequency. The frequency-dependent absorption opacity is denoted by k^ . In general, the functional dependence of a 
given quantity is suppressed unless it is necessary to distinguish between functions with different arguments (e.g., equation 1111 '). 

We consider both the effects of absorption and scattering, but treat them in an approximate manner in order to ensure algebraic tractabil- 
ity. For simplicity, we include scattering only in the shortwave via the param eter £, whic h is the ratio of the absorption to the total opacity, 
while being aware that longwave scattering is a non -negligible eff ect ( de Kok et al.ll201in rl Thus, we have ^ = and 1 in purely scattering 
and absorbing situations, respectively. For example. iGuillod J2010l) assumes ^ = 1. It is important to note that decreasing ^ has the effect of 
increasing the number of shortwave scatterers present, while keeping the number of shortwave absorbers constant. The Bond albedo A and 
^ are not independent parameters. Rather, as a generic statement, we expect: when ^ = 0, .4 = 1; conversely, when ^ = 1, ^ = 0. These 
properties serve as boundary conditions to the functional form ^ = ^(.4), which we will elucidate later, in the context of the coUimated beam 
approximation (Appendix IaTi. in equation dlVb and Figure (TJ. 

The shortwave optical depth is 



TS 



— - dm 



Ksm 



(1) 



/o <, C ' 

where m is the column mass per unit area and ks is the (constant) shortwave absorption opacity which we will describe shortly. In a 
hydrostatically balanced atmosphere, the shortwave photosphere resides at 

p _ ^»? 

with g being the surface gravity of the hot Jupiter. The preceding expression follows from demandin g that ts = zo.ln the solution to the 
classical Milne's problem for a self-radiating (stellar) atmosphere (see Section 3.3 of lMihalaall978h . one obtains vj = 2/3 ~ 0.67. For 
irradiated atmospheres, the value of vj needs to be estimated numerically (see Appendix |BJ; surprisingly, its value is not very different: 



(2) 



The formalism remains algebraically tractable even if the longwave analogue of § is included as an additional parameter. 



© 2011 RAS, MNRAS 000,[T]tT8 



Clouds and hazes in hot Jovian atmospheres 3 

■uj ~ 0.63. The photon deposition depth — the location where the energy of the stellar photons are deposited — resides slightly deeper 
than the shortwave photosphere: at \/S,ts — zu instead of rs — vj (see equation 1151 ). Thus, the pressure level at which shortwave photon 
deposition occurs is 



Pd = 



63 mbar y^ 



(3) 



Ks V4 ■ \ 0.63 10 m s-2/ V 0.01 cm2g- 

As the effect of scattering becomes stronger (decreasing 5), the altitude at which the stellar photons are absorbed becomes higher if shortwave 
absorption stays constant. In the extreme limit of ^ = 0, no starlight gets absorbed and the atmosphere has an albedo of unity — its 
temperature-pressure profile is described only by its internal heat flux. Since ^ ^ 1, we have Pd ^ ^s. When shortwave scattering is absent, 
we obtain Ps = Po- 

For the rest of the paper, we refer only to absorption (and not scattering) opacities when we use the term "opacity" (unless otherwise 
specified) and assign symbols to them. The shortwave opacity, 

_ Jg K^ J^ df 



/g J„ du 



(4) 



is assumed to be co nstant with the integration being taken over the range of shortwave frequencies. Equation (HJl defines the absorption mean 
opacity (see §3.2 of iMihalaJI 197a) . which is defined to guarantee the correct total amount of (shortwave) energy absorption. By contrast, the 
longwave opacity, 

/j^ KvJv dp J-^ HvBv dv 



KL 



j^ Jv dv J^ B„ du 



(5) 



is closer to bein 



; a Planck mean opacity defined to guarantee radiative equilibrium in the optically thin portions of the atmosphere (§3.2 



of iMihalasll 19781) . The Planck/blackbody function is denoted by B^. The longwave opacity is in general a function of m, the cosine of the 
latitude n = cos 6 and the longitude (j>, i.e., kl = i^Lim, /i, (p). For simplicity, we consider kl = i^Lim) only, while acknowledging the 
possibility that the strong temperature gradients inherent in hot Jupiter atmospheres may result in latitudinally and longitudinally varying 
longwave opacities. 

Denoting the specific intensity by /„ (with units of erg cm^'^ s~^ Hz^^ sr~^), the radiative transfer equation is (e.g., pg. 35 of 

1978h 



Mihalas 



Ml 



ulu Kii/li^ 



t^i/Di/ 



OJ^ 



(6) 



dm e ■'"'""' e 

Equation ^ implicitly assumes either isotropic scattering or the collimated beam approximation (see Appendix |AI where the for ward and 
backward scattering probabilities are the same. The mom ents of the radiative transfer equation are described by a pair of equations JMihalas 
19781: iHubenv. Burrows & SudarskvlboOslJGuillotlboiol) , 

^^•^ — f T R ^ 

(7) 



dm 



l^i^Hu 



dm ^ 

Since there are three unknowns {Jv,Hi,,Kv) and only two equations, closure relations known as the Eddington approximations are required 
in order to obtain the solutions. In the longwave, the first and second Eddington coefficients are respectively defined as 



£i 



(8) 



Jl 3 Jl 2^ 

and are assumed to be constant — with their values chosen to be consistent with lGuilloll J2010l) — to facilitate algebraic amenability. We 
note that if the two-stream a pproximation is made — wh ich is typic ally the simplifying assumption adopted in three-dimensional simulations 
with radiative transfer (e.g. Jshowman et alj |2009: Hens. Frierson & PhillipDS 2011) — then we instead have, e.g., £^2 = l/%/3 ~ 0.58 (see 
Appendix IaI). In the shortwave, the first Eddington coefficient is equated to the square of the cosine of the latitude ( lGuilloll20ld) . 

Js' 



^^■ 



(9) 



We note that equation Q is a consequence of the collimated beam approximation which we make only in t he shortwave (see AppendixlAt 



_ , _ . PP e 

Let the heat transported by atmospheric circulation, as witnessed in three-dimensional simulations (e.g.. lShowman et al.l2009luHeng. Menou & Phillippsl 
Q{m, fi, (f>), which has units of erg s~^ g~^ and is thus technically a specific luminosity of heat. It can be related to the mo- 



20111) . be O 

ments of the specific intensity by integrating the right-hand side of the first expression in equation Q over all frequencies, 

kl(Jl-B) + ksJs = Q, (10) 

where we have B = j^^B^ dv = osbT'^'/tt, ctsb denotes the Stefan-Boltzmann constant and we have assumed that the star and the hot 
Jupiter radiate negligibly in the longwave and shortwave, respectively. Equation l llOt is equivalent to the energy equation used in general 
circulation models if one specializes to a static or time-independent situation. By further integrating over column mass, it follows that 

H = H^~Q{m,oo), (11) 
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Figure 2. Examples of the first moment of the shortwave intensity for ^ = 1 and Ti^j = 2025 K. 



where Hoo is the value of H evaluated when m -^ oo and 

Q (7711,7712) = I Q (m' , /i, (j)) dm'. 



(12) 



Note that equation l lllb is completely general: it requires no assumptions on the functional forms of either ks or kl- For convenience, the fi- 
and (^-dependences of Q have been suppressed. 



2.1 Shortwave 



In the shortwave, equation lO reduces to 



dHs 
dm 
dKs 
9m 



= KS Js, 
usHs 



(13) 



C 



By using ^^ — Ks/Js, one may rewrite the preceding pair of equations. 



& 



'Js^ ( ^s y 



dm? \^.y/lJ 

3rr / . \ 2 



Js + 



d- 



1 d{ns/0 
1 dns 



KsJs dm' + Hso 



dm^ X^yJXJ y? dm 



KsHs 



dm' + ^i^Jso 



(14) 



/o C 

Upon inspection of equation J14t . it is evident that obtaining analytical solutions for Jg and Hs is a challenging task unless the assumptions 
of constant ks and ^ are made. With such assumptions and the boundary conditions that Js = Hs = as m — >■ cxd, we obtain 

Ksm\ f y^Ts'' 



Js = Jso exp 



^s = Hso exp 



(15) 



where Jsq and Hsg are the values of Jg and Hs, respectively, evaluated at m, = and whence 

Hs = -f^\/iJs => Hso = -fJ'VlJso- 



(16) 



One may regard Kg as a shortwave opacity which collectively describes the absorbing effects of the "gas" and the "haze" or "cloud", since it is 
assumed to be constant. Thus, the effects of absorption and scattering in the shortwave are described by the parameters ks and ^, respectively. 
Combining equation J16b with the coUimated beam approximation already employed in the formalism (see Appendix |A](, the Bond 
albedo may now be related to the scattering parameter (Figure[T](, 



A- 



(17) 



1 + %/e' 

which in turn allows us to make a physical connection between the strength of shortwave scattering and the temperatures deep in the 
atmosphere (Toe; equation 1441 ). Note that we only consider forward and backward scattering in the coUimated beam approximation (see 
Appendix [At. 

As an illustration, we show examples of Hs in Figure [2l Two effects are apparent: the stellar photons are absorbed higher in the 
atmosphere as scattering becomes more important; and stronger scattering dilutes the incident stellar irradiation via the enhancement of the 
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Bond albedo. Both the shortwave photosphere (equation (2)) and photon deposition depth (equation (51) shift to higher altitudes as scattering 
becomes a stronger effect (decreasing ^), but the relative separation between them grows as oc l/\/C- 



2.2 Longwave 

In the longwave, equation Q reduces to 



dm 

dm 



= kl{Jl — B) = Q — KsJs, 



(18) 



where the second equality in the first expression follows from the assumption of energy conservation as described in equation JlOt . Integrating 
over the column mass (per unit area) for the first expression in equation l llSt . we obtain 



-ffL = -ffLo + Hso 



exp 



KgTTl 



+ Q(0,m). 



(19) 



In deriving the preceding expression, we have made use of equation ( I16t which assumes a constant shortwave opacity. However, equation 
( 119b makes no assumption on the functional form of kl. The quantity H-l^ is the value of Hi^ evaluated at m = 0. Using equation dllb . it 
follows that 

Ksm\ 

Combining equation l[8), the second equation in l llSb and equation i20b yields 



Hl ~ Ha 



Hso exp 



• (m, oo) . 



(20) 



t-1 Jo 



-ffoo - Hs(, exp 



Ksm, 



:\ 



' (m,oo) 



dm , 



where Jlq is the value of Jl evaluated at m = and may be eliminated in favour of other quantities using equations ([8} and MQ\ . 



Iffco - Hso - Q (0, oo; 



(21) 



(22) 



Our goal is to obtain an expression in terms of only B, Hoo and Hsq ■ To this end, equation l l21t may be rewritten using equation l llOt . 



B = Ha. 



1 



kl dm 



Hs 



1 1 

£2 £1 JO 



kl exp 



1 

^£2 £1 JO 

where we have collected all of the terms involving Q into the following quantity, 

1 



Ksm'\ I 
-^r dm 



+ 



«:s 



M\/?'^L 



exp 






+ Q, 



K-L tl Jo 



KL 



' (m',00) dm,' — 



'(0,00) 
£2 



(23) 



(24) 



We have defined Q < because it is the heat flux — as a function of depth, latitude and longitude, i.e., Q — Q{m, /i, (jj) — transported by 
horizontal winds from the day to the night side of a hot Jovian atmosphere ( Showman & Guillot 2002) . In the li mit of a constan t longwave 
opacity (i.e., dKi^/dm — 0) and no shortwave scattering (^ — 1), equation ( 123b reduces to equation (41) of iGuilloll ( 120100 when two 
typographical errors in the latter are corrected: one should obtain — Q/kl instead of — Q and the terms involving Q{m' , 00) are missing the 
£^^ coefficient. 



2.3 Global-mean temperature-pressure profile 

As realized bv lGuilloll ( l2010h . the pair of first moments in equation J23b have clear physical interpretations: 

fSBrint rr M'^SsTirr 



Ha 



Hso 



4n ' " 47r 

The quantity Tint is the blackbody-equivalent temperature associated with the internal heat flux, while the irradiation temperature is 



-'-ITT ^ * 1 

a 



1/2 



-A) 



1/4 



1900 K 



f r* 



V6000K 



IL/o 
0.1 



1/2 



2%/C 



i^ + VO 



1/4 



(25) 



(26) 



where T* is the stellar effective temperature, _R* is the stellar radius and a is the spatial separation between the hot Jupit er and the star . The 
albedo integrated over all shortwave frequencies is denoted by A and may be regarded as the Bond albedo (see §3.4 of lSeageauOld for a 
discussion of albedos). Note that the equilibrium temperature, in the absence of day-night heat redistribution, is Tcq — Tirr/v^- The quantity 
Hso is negative because the incoming stellar irradiation travels downwards into the atmosphere. For most values of the scattering parameter 
(5 > 0.1), the coefficient (1 — .4)^'* = [2-^/(1 + \/C)]^^* is typically close to being unity (Figure[T](, which implies that the temperature 
at depth Too is insensitive to changes in the Bond albedo (cf. through Tcq in equation 1441 ) unless A > 0.5. 
Substituting equation l l25b into equation l l23b . we obtain 



r" 



T. 



1 



^ ( ^ + ^ / KL dm' ) + ^ 



/y-i4 






exp 



Ksm\ , _/^ I _M_ 
mVC/ £2 £1 JO 



KL exp 



Ksm 



n 



t^V^J 



dm 



+ -^. (27) 

crsB 
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In the limit of a constant longwave opacity, the preceding expression reduces to equation (43) of lGuillotI d201(]|) when a typ ographical in the 
latter is corrected for: instead of TvQ/asB, one should have — 7rQ/(7SB«:L. The reduction to equation (43) of lGuilloll OOICI) also requires the 
use of the identity (see Appendix [Q. 

d 



dm 



■Q (m , oo) 



•{m',^,(f>) 



(28) 



X= — 

27r 



towards evaluating the second term in equation l l24t via integration by parts. 
For an arbitrary function X, its global-mean counterpart is 

-| rllT rl 

/ / Xdiidij), (29) 

/o Jo 

with the assumption that there is latitudinal symmetry about the equator (6 — 0°). For the terms involving the internal heat flux (o-gBTifj^), 

we assume azimuthal symmetry, i.e., the day and night sides are the same. For the terms involving the irradiation temperature (Tirr), we 

are strictly speaking performing only a hemispherical averaging — the integration over is only non-zero for vr radians over the day side, 

consistent with the theoretical expectation that hot Jupiters are tidally locked (at least in the case of circular orbits). If one wishes to compute 

the temperature-pressur e profile only on the day-side of a hot Jovian atmosphere, then the integration over needs only to be performed 

from to 7r|j Following I GuillotI feoioh . we assert that 

-1 r2n rl 

— Qdfid(l> = 0. (30) 



The three-dimensional simulations of atmospheric circulation of iHeng. Frierson & Phillippj ( 120111) demonstrate that the approximation in 
equation i30b breaks down at P < 10 bar, where the simulated temperature-pressure profile (in dynamical-radiative equilibrium) is cooler 
by ~ 50 K, at least for their models of HD 209458b, due to the conversion of heat into mechanical energy via the creation of ~ 1 km 
s^^ horizontal winds. Since our goal is to derive temperature-pressure profiles in radiative equilibrium and not full dynamical-radiative 
equilibr i um, w e retain this approximation. One may regard our approach as being analogous to the "no redistribution" models presented in 
HansenI hoOSh . Thus, we get 



T* 



4 



1 1 



HLh dm ) + 



-^ ir 






Ksm 
where we have explicitly stated the arguments of the exponential integrals and defined 



+ 



£i 



KL Es 



Ksm 



dm 



7^ 



KL 



'-, T — T (m) = Ki^m. 



(31) 



(32) 



We call the quantity T the "global-mean temperature", but strictly speaking it is the 1/4-root of the global-mean flux since the averaging is 
taken over T* and not T. The exponential integral of the j-th order is described by the following expressions (lArfken & Weber 1995) : 



Ej (x) = 
Ej+i (x) 



y ■' exp i-xy) dy, 



(33) 



[exp( 



xEj (x) 



If we set kl to be constant, us e S, = 1 and ad opt the numerical values of the Eddington coefficients as stated in equation {8](, then e quation 



3\\ re duces to equation (49) of lGuilloll ( 1201 Of) . Note that none of the typographical errors previously described affect equation (49) of lGuillot 
20101), because of the assertion in equation ( 1301 ). 



In summary, equation ( 131b is the general expression for the global-mean temperature, as a function of depth, in a hot Jovian atmosphere, 
assuming constant shortwave opacity but allowing for an arbitrary functional form for the longwave opacity. In the following section, we will 
examine various forms of this equation with kl = i^him,) specified. 



3 APPLICATIONS OF THE FORMALISM 

To proceed further, we first need to relate the column mass (per unit area) m to the pressure P. If we specify z to be the vertical 
coordinate as measured from the top of the atmosphere (at z — 0) downwards, then hydrostatic balance demands 

dP 



dz 



^ pg, 



(34) 



where p is the mass density of the atmospheric fluid and g ~ 10"^ cm s ^ . By assuming that g remains constant over the extent of the (thin) 
atmosphere and defining dm = pdz, we obtain 



P — mg. 



(35) 



In other words, the term in equation )3U involving T^^^ has a coefficient of 1/8 and not 1/4. 
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Figure 3. Temperature-pressure profiles with 5 = 1 and various values of e, which parametrizes the effect of collision-induced absorption in the longwave. 



For the purpose of comparison to three-dimensional simulations, which require the specification of a bottom (with a pressure Po) for the 
computational domain, we define the column mass per unit area at the bottom to be 



mo = Pol 



(36) 



Therefore, we have P/ Po = m/mo- We adopt Po — 100 bar throughout this study. 

Upon specifying a functional form for kl = kl(7Ti), the simplification of equation ( I31t only involves the evaluation of the double 
integral. 



-^y.-v. -M-(-^i-- 



For example, when kl is constant, we have 



1 = 






l-3i..(I^ 



(37) 



(38) 



3.1 Collision-induced absorption 

As previously noted in lGuilloti ( 120 101) and iHeng. Frierson & Phillippsl d201lh . the assumptio n of a constant l ongwave opacity breaks 
down at high press ures because of the effect of collision -ind uced absorption JHer zberg 1952; Pi errehumbertll201Q) rl We then have kl oc 
P cc m. Following lFrierson, Held & Zurita-Goton ( 120060 and lHeng. Frierson & Phillippsl ( 120110 , we adopt the following functional form for 
the longwave optical depth. 



TL (P) = TO 



W.*^'-<K 



(39) 



such that tl ~ toP/Po and tl = ero at the top and bottom of the model atmosphere, respectively. The dimensionless quantity e is the 
correction factor to the longwave optical depth due to the effect of collision-induced absorption, which becomes important when P > 
Po/(e — !)• The functional form for the longwave opacity is 



kl (m) = Ko 1 + 2 (e - 1) 

which guarantees that equation ( |39t is recovered when we evaluate 

TL (m) 



mo 



Ko, m < mo, 

Ko (2e - 1) , m = mo, 



kl [in 



') dm' 



(40) 



(41) 



^ Note that colhsion-induced absorption and pressure broadening are different effects: the former produces continuum absorption, while the latter is associated 
with the broadening of absorption lines. 
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Figure 4. Hypothetical cloud/haze decks parametrized using equation j46t . Shown are three examples with different values of the deck thickness parameter. 



and use m/mo = P/Pq. The normalization for the longwave optical depth in the absence of collision-induced absorption is tq 
follows that equation i3U reduces to 



y4^37^ 



+ TL + 



m, 



+ 



STeljC 



^o7o 



exp 



M 
370 

'7! 



1 + cxp 



3-fT 



£4 



-y-7- 

7! 
7! 






27 /7T 



1- 



Ko'TT-o. It 



(42) 



where the functional dependences of the exponential integrals on 7t/v^ have been explicitly stated and 



70 = — 






(43) 



In the absence of clouds/hazes, atmospheres with and without temperature inversions are characte rized by 70 > 1 and 70 < 1, respectively. 
When e = 1 and ^ = 1, we obtain tl = r, 7 = 70 and equation ( I42t reduces to equation (49) of lGuillotI OOlOl) . The temperature at depth 
(Tint = K, T > 1) is 



Tx 



:Te< 



1 I ^^ I 3^(^-1) 

2 ^ 270 



-,1/4 



470 rso 



(44) 



The quantity tsq is the value of the shortwave optical dept h, ts = Ksm, evaluated at m = rwp. Figure [5] shows examples of temperature- 
pressure profiles, with ^ = 1, adopting the numbers used in lHeng. Frierson & Phillippj ( 1201 Ih for their models of HD 209458b: Tint = K, 
Teq = 1432 K, Ks = 0.006 cm^ g~^, acq = 0.01 cm'^ g~^ (70 — 0.6) and tsq = 1401. It is apparent that larger values ofe correspond to 
higher temperatures at depth, an effect which is similar to line blanketing in stellar atmospheres. For e = 1000, the temperatures at P ~ 0.01 
bar appear cooler than the other cases because the enhanced opacity contribution from collision-induced absorption becomes noticeable at 
these pressures, similar to the cooling effect due to the presence of a cloud/haze layer (see i]4.1.1t. Elucida ting the chemistry/physics which 
determines the exact value of e is beyond the scope of the present study (see §4.4.8 of IPierrehumberlll 20 1 d for a review of collision-induced 
absorption). Rather, we demonstrate that the effect of collision-induced absorption can be approximately described using a a single parameter. 
Note that the value of Too obtained, for an assumed value of e, depends on the value of Pq adopted — in other words, one needs to know 
the enhancement of the longwave optical depth, due to collision-induced absorption, at a reference pressure level (e.g., To)- Another way of 
interpreting the results in Figure[3]is that if e < 100 at Pq ~ 100 bar, then collision-induced absorption is a min or effect (i.e., temperature 
differe nces < 100 K). As furtherexamples, we re-compute the values of Too, for^ = 1. used in the simulations o fiHeng. Frierson & Phillippj 
i 201lh : for Po = 220 bar and setting e = 10, 100, 1000 and 2000, we h ave Too j ^a 154 1, 1564, 1749 and 1903 K, respectively. 

Finally, we not e that the (simpler) general i zation of Section 2 of lGuilloll ( 120100 to consider equations ( |39b and J40b was previously 
presented in §4.3 of Heng. Frierson & Phillippj ( 120 111) . Also, the approximations of constant ks and kl oc P were previously adopted b' 
Liu & Schneiden ( 120100 for simulating the atmospheres of the gas giants in our Solar System. Furthermore. lHeng. Frierson & Phillippj j201ll 
used the functional form in equation {i9\ to study the atmos pheres of hot Jup i ters an d were able to produce zonal-mean flow quantities (e.g., 
wind, temperature) which are similar to those published bv lShowman et al.l 120091), who used full opacity tables in their multi-wavelength 
calculations. 
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3.2 Uniform cloud/haze layer with absorption and scattering 



In the simple case of a uniform cloud/haze layer, one can include the effects of both shortwave scattering (^ < 1) and an extra, longwave 
opacity contribution Kcq ■ Equation i31t then becomes 
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(45) 



where we have 7c = ks/kcq • Equation l l45t is essentially an analytical expression short of having to evaluate the exponential integrals. 

In the longwave, we have made a distinction between the extra cloud/haze opacity Kcq and the opacity kl describing the rest of the 
atmosphere (the "gas") because we assert the former to be constant while the latter is subjected to collision-induced absorption. One may 
also choose to set Kcq ~ cm^ g^^ and simply assimilate the longwave absorption properties of the haze/cloud into kl. In the shortwave, 
the opacity contribution due to the cloud/haze can likewise be added to that of the gas because the shortwave opacities are assumed to be 
constant. 

3.3 Purely absorbing cloud or haze deck 

^The formation and dis persal of clouds or hazes is a fundame ntal, unsolved problem in the study of terrestrial climate science JPierrehumberd 

2010l) . brown dwarfs (e.g.. lBurrows. Heng & NamDaisanj|201in and exoplanetary atmospheres. To keep our formalism semi-analytical, we 
allow for the presence of a non-uniform cloud deck via the (ad hoc) function. 



Kc — ^CQ 6Xp 



-Ac 



-^ 



(46) 



Thus, our cloud deck is described by three free parameters: the opacity normalization Kc,,, the deck thickness parameter Ac and the location 
of the deck P^. Figure |4] shows three examples of hypothetical cloud decks located at Pc = 0.1 bar and with different deck thicknesses — 
it is apparent that the larger the value of Ac, the thinner the cloud/haze deck. In the limiting case of Ac = 0, equation l |46t describes an 
extra opacity contribution which is constant throughout the atmosphere. In the interest of algebraic amenability, we consider the cloud/haze 
deck to contribute an extra source of opacity only in the longwave; we set ^ = 1 since a constant scattering opacity is inconsistent with 
an extra, non-uniform absorption opacity. While being a simpl e starting point for non-uniform decks, real clouds or hazes probably do not 
exhibit this type of simple behaviour. (See lPierrehumberlluOlOl for a review of the physics of clouds an d their effects on the atmosphere.) For 
example, clouds in brown dwarfs have been observationally inferred to be patchy JArtigau et al.ll200 9"). Nevertheless, equation J46t provides 
us with a lucid vocabulary for describing cloud/haze decks: dense (higher Kcq) versus tenuous (lower Kcq), thick (higher Ac) versus thin 
(lower Ac), high (lower Pc) versus low (higher Pc). The main motivation in this sub-section is to isolate the effect of an extra, non-uniform, 
infrared/longwave source of opacity. 

The additional longwave optical depth due to the cloud/haze deck is 



Kc dm 
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.1/2 



m 
nic 



where we have defined rric 



Pc/g and the error function is ( JArfken & Webeill 19950 . 

E{x)= —= / exp (-y^) dy. 
V^ Jo 



Generalizing equation ( |39t , the longwave optical depth becomes 



TL = To 
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It follows that the global-mean temperature-pressure profile becomes 
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with the integral J being described by 
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(51) 



where the arguments of the exponential integrals and error functions have again been explicitly written out. The integral J does not have a 
general analytical solution, needs to be evaluated numerically and accounts for the warming effect of the clouds in the lower atmosphere. 
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Figure 5. Global-mean temperature-pressure profiles in the limiting case of constant cloud opacity (Kc = Kcq ) and 5 = 1. Left panel: Kg = 0.006 cm^ g ^ 
(70 = 0.6). Right panel: Kg =0.2 cm^ g~^ (70 = 20). The equilibrium temperature of the hot Jupiter considered here is Tcq = 1432 K. 



3.4 Longwave opacity with power-law term 

For completeness, we note that if the longwave optical depth and opacity have functional forms consisting of linear and power-law 
terms, 

n + l" 



TL (P) = TO 



kl (m) = Ko 



P_ 



(^-1) 



Po 



l + (n+l)(e-l) 



/ m 



where n > is dimensionless, then the simplification of equation l l31t requires the numerical evaluation of the integral, 
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4 SPECIFIC EXAMPLES OF TEMPERATURE-PRESSURE PROFILES 

All of the results in this section assume Tint = K. 

4.1 Absorption only 

All of the results in this sub-section assume ^ = 1. 



4.1.1 Profiles with constant cloud/haze opacity 

Using equation l l45b . Figure [S] shows examples of temperature-pressure profiles for Kc = Kcq snd adopting parameters similar to those 
for the hot Jupiter HD 209458b: T^q = 1432 K, g = 9.42 m s'^, kq = 0.01 cm^ g-^ For illustration, we have set yi = and e = 100 
— the former follows immediately from ^ = 1 and equation J17l (. while the latter implies that collision-induced absorption is the dominant 
source of the longwave optical depth at P > 2 bar. Higher values of e simply translate into higher values of the temperatures at depth Too 
(Figure[3}. 

The left and right panels of Figure[5]show the cases of 70 = 0.6 and 20, respectively. The effect of the cloud/haze, present throughout 
the atmosphere, is to cool and warm the upper and lower atmospheres, respectively. As the cloud opacity (kcq) is increased, the longwave 
photosphere located at Pl ~ g/{i^o + Kcq) shifts to higher altitudes. For 70 = 0.6 (left panel), the lower atmosphere (P > 0.1 bar) has a 
temperature which is generally higher than the equilibrium temperature of Tcq = 1432 K — the analogue of the greenhouse effect for hot 
Jupiters. For 70 = 20, the effect of the cloud/haze is still to cool and warm the upper and lower atmospheres, respectively, but the lower 
atmosphere now has a temperature which is generally lower than the equilibrium temperature — the analogue of the anti-greenhouse effect 
for hot Jupiters. (See §4.3.5 of |Pierreh umber][ 201CI for a discussion o f the greenhouse and an ti-gr eenhouse effec ts.) 

Previously, it was pointed out bv lHubenv. Burrows & Sudarskvl <2003r ). lHansea ( 120080 and lGuilloll ( 120 10 1) that temperature inversions 
occur in an atmosphere when 70 > 1. When a uniform cloud layer is present, the condition for a temperature inversion being present becomes 



70 > maxjo, (l-7c ^) ^| • 
In other words, if 7c > 1, then whether 70 > 1 becomes irrelevant. When clouds are absent, we have 7^ 



(54) 
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Figure 6. Global-mean temperature-pressure profiles with purely absorbing Gaussian cloud/haze decks (assuming ^ = 1), which illustrate their dual effects of 
cooling and warming the upper and lower atmosphere, respectively. Shown are decks with three different thickness and two different opacity normalizations. 
The left and right columns are for Pc = 0.1 bar and 0.01 bar, respectively. The top, middle and bottom rows are for 70 = 0.6, 2 and 20, respectively. In all of 
the panels, the cloud-free case is shown for comparison. 

4.1.2 Profiles with Gaussian cloud/haze decks 

In Figure|6l we evaluate equation J50b numerically to obtain the temperature-pressure profiles for models with various realizations of the 
idealized Gaussian cloud deck. We retain the values of Tcq, g, kq, ng, A and e previously described in i]4. 1.11 Unlike in the case of a uniform 
cloud layer, the cooling/heating effect of the cloud deck depends on its location with respect to the shortwave photosphere as described in 
equation ^. For Kg = 0.006, 0.02 and 0.2 cm^ g~^, we have Pg « 0.1 bar, 30 mbar and 3 mbar, respectively. 

In the top row of Figure |6l we adopt us = 0.006 cm^ g~^ such that 70 = 0.6. In the top left panel, the cloud deck is essentially 
coincident with the shortwave photosphere. Similar to the case of a uniform cloud layer, the cloud deck cools and warms the upper and 
lower atmosphere, respectively. When the cloud deck is moved upwards to Pc — 0.01 bar (top right panel of Figure [6](, above the shortwave 
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Figure 7. Global-mean temperature-pressure profiles with uniform, scattering cloud/haze layers, which illustrate their dual effects of cooling and warming the 
lower and upper atmosphere, respectively. The left and right panels are for Kcq = and 0.01 cm'^ g~^, respectively. In both panels, the no scattering case 
(5 = 1) is shown for comparison. 



photosphere, the magnitude of the cooling effect remains somewhat the same but the warming effect becomes less pronounced. In addition, 
the thinner cloud decks (e.g., Ac = 10) imprint their Gaussian shapes onto the temperature-pressure profiles, causing a small temperature 
inversion. 

When the cloud deck is located far below the shortwave photosphere, as shown in the bottom left panel of Figure [6] where we adopt 
Ks = 0.2 cm^ g~^ such that 70 = 20, it still cools the upper atmosphere but has no effect on the lower atmosphere (P > 0.02 bar). 
Furthermore, only when the cloud deck is both thick (lower Ac) and dense (higher Kc,,) are the temperatures noticeably cooler compared to 
the cloud-free case. When the cloud deck is shifted higher in altitude to Pc — 0.01 bar, the warming of the lower atmosphere now becomes 
noticeable, but is still more subdued than in the case when Pc < Ps- 

The middle row of Figure[6]possesses the most physical richness compared to the other two rows. The shortwave photosphere is now 
located at Ps ~ 30 mbar, which corresponds to 70 = 2. In the absence of a cloud deck, one would expect the temperature-pressure profiles 
to always possess temperature inversions. For a Gaussian cloud deck, even the generalized condition in equation i54t is insufficient to predict 
if a temperature inversion will occur. When the cloud deck is tenuous (kcq = 0.01 cm^ g~^), thin (Ac — 3, 10) and low-lying (Pc — 0.1 
bar), temperature inversions occur, consistent with equation l l54t . However, if the cloud deck is made to be denser (Kcq = 0.1 cm^ g~^) 
it introduces a component to the profile where the temperature decreases with increasing altitude. If it now sits higher in the atmosphere 
(Pc — 0.01 bar), then there are two inverted components in the profile along with a non-inverted component. By contrast, a low-lying 
(Pc — 0.1 bar), thick (Ac = 1) and dense (kcq = 0.1 cm^ g~^) cloud deck "reverses" the 70 > 1 condition and produces a "normal" 
temperature-pressure profile which behaves as if 70 < 1. All of these features illustrate the complexity introduced by a simple, Gaussian 
cloud deck and demonstrate that there is no straightforward way to generalize the condition in equation i54t . 



4.2 Inclusion of shortwave scattering 

We now allow for ^ 7^ 1 and a non-zero Bond albedo (via equation 1171 ) and investigate the effects of shortwave scattering by a uniform 
cloud/haze layer. Setting Kcq = cm^ g~ ^ , the left panel of Figure|7]demonstrates that shortwave scattering has the opposite ejfectfrom that 
of an extra longwave opacity: it heats and cools the upper and lower atmosphere, respectively. This is further illustrated by comparing the 
left and right panels of FigurelT] where one sees that "turning on" the cloud opacity Kcq partially negates the effects of scattering. In general, 
scattering either produces or strengthens a temperature inversion as the shortwave photosphere is shifted to higher altitudes. 

So far, we have used the terms "cloud" and "haze" synonymously. We now attempt to better define the two terms. By "cloud", we refer 
to a layer of enhanced optical depth caused by certain chemical species condensing out of their gaseous forms. In practice, these cloud decks 
will be situated roughly at the pressure level where their condensation curves intersect the temperature-pressure profile of the atmosphere. 
Thus, clouds may be regarded as having a thermody namic origin. Such an approach of "p ainting on" a cloud deck has been used to construct 
spectral-evolutionary models of brown dwarfs (e.g.. lBurTows. Heng & Nampaisarnll201 11) . By contrast, a "haze" layer may refer to an extra 
source of (absorption and scattering) opacity which is either uniform throughout the atmosphere or is located at a pressure level which does 
not depend on a particular condensation curve. Thus, hazes may be regarded as having a non-thermodynamic (e.g., photochemical) origin. 
The simplicity of our models allows us to explore the effects of clouds and hazes on the temperature-pressure profile without having to 
perform a complicated coupling to a thermo-chemical model; the latter is required not only for self-consistency but also if predictions for 
chemical species and particle sizes are desired. 
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5 DISCUSSION 

5.1 Application of models to observations: HD 189733b 

We now make a tentative attempt to compare our models of the global-mean temperature-pressure profile to data points inferred from 
the published observations of the hot Jupiter HD 189733b. We focus on HD 189733b as a case study because of the detection of haze in its 
atmosphere (Pont et al. 2008; Sing et al. 2009, 201 1). We make use of the transit spectroscopy, secondary eclipse and photometr ic phase curve 



observ ations of HD 189733b from the Hubble and Spitzer Space Telescopes. In comparing one- to three-dimensional models. iFortnev et al 



1 2010h demonstrated that the temperature-pressure profiles on the exoplanetary limb are good approximations to the true global-mean profile, 
making it natural to apply our formalism to temperatures derived from transit data. 

We extract the hemisp herically-averaged brightness temperatures from the phase curves measured using Spitzer at both 8 and 24 /im 
(JKnutson et al. 2007, 2009). For the temperatures at the limb during transit, we use the average value, (Tmax + Tmin)/'^, as given in Table 3 
oflKnutson et al. (2009) for both the 8 and 24 /im temperatures, fin ding temperatures of 1135 =fc52 K and 1102zfc67 K, r espectively. We also 
use the 3.6, 4.5, 5.8 and 16 ixm secondary e clipse measurements of ICharbonneau et al.l ( 120080 and lDeming et al.l j2006h with the brightness 



temperatures given in lCowan & Agoll 1 1201 ih . As there is no phase curve information at these wavelengths, we applied a correction of one-half 



of the average day-night temperature contrast as measured at 8 and 24 ^m (121±35 K) to estimate limb temperatures, finding temperatures 
of 1518 ± 49 K, 1197 ± 57 K, 1247 ± 77 K and 1217 ± 62 K at 3.6, 4.5, 5.8 and 16 /x m, respectively . The a ssumed pressures for these data 
points are model dependent and based upon the normalized contribution functions from Knutson et al.l (120091) . As such, the absolute pressure 
scale is uncertain and ultimately tied to the assumed (solar) model abundances and composition, with H2 O typically bein g the dominant 
opacity source in the near infrared. Water has been identified in absorption in the Spitzer emission spectra of lGrillmair et alj ( 120081) . 

We augment the Spitzer data with the temperatures associated with Rayleigh-scattering haze as detected and confirmed by Hubble 
Space Telescope (HS T) transmiss i on sp ectra ( font et al . 2008; Sing et al. 2009, 2011). We use the temperature derived from the HST ACS 
and NICMOS data in Ising et alJ ( 120091) of 1280±110 K. We tie the pressure of the optical haze to that of the Spitzer emission spectra 
measurements using the differences in altitude between the optical and near-infrared transit spectra as well as the difference from the transit 
ge ometry. Denoting the radius of the hot Jupiter by Rp, the 0.75 ^m optical transit radius is 7 x 10~*-Rp/-R* above the 8 /im transit radius 
of lAgol et al.l ( I2010I) . which is approximately 1.7 scale heights at 1300 K or a difference in pressure of a factor of 5.5. In addition, the 
slant-to-normal difference of the transit geometry decrease s the pressure of the transit measurements by a factor of ^2-KRp/ H ~ 50 at a 
given wavelength, with H being the pressure scale height JBurrows et al . 2001; Fortney 2005; Hansen 2008). Thus, the haze in the optical 
at 0.75 fji,m is a factor of about 275 lower in pressure than the 8 /^m point from the emission spectrum, placing th e haze at about the 0.7 
mbar pressure level. This estimated pressure value is in agreement with that of iLecavelier des Etangs et alj (|2008J), who assumed a haze 
composition of MgSiOa. 

While P < 0. 1 mbar data points do exist (Huitson, Sing et al., in preparation), it is important to note that o ur formalism is invalid at these 
low p ressures because it does not include the physics of the thermosphere, including hydrodynamic escape ('e.g.. lMurTav-Clav. Chiang & Murray 
20091) . It is likely that our assumptions of hydrostatic balance and radiative equilibrium break down significantly at these altitudes. Therefore, 
we terminate the model comparisons for P ^ 0.1 mbar. 

T he initial task is to cut do wn on the number of parameters used in our models. We estimate that Tirr ~ 1697 K (1—^)^'* dBouchv et alj 
20050 and g « 21.88 m s~^ dTorres. Winn & Ho lman 2008). Since the data points are for P < 1 bar, they set no constraint on the 
temperatures at depth and thus are unaffected by the value of Tint adopted — for simplicity, we set Tint = K. Again for simplicity, we 
set KcQ = cm^ g~^, namely that any extra longwave absorption by the cloud/haze present may be assimilated into the longwave opacity 
function kl. Therefore, we are left with a 4 -parameter mode l to be compared to 7 data points. The assumption of a uniform scattering 
opacity is not unreasonable given the fact that lSing et al.l J201lh detect a transmission spectrum of HD 189733b consistent with a haze layer, 
producing Rayleigh scattering, at ~ 1 mbar, spanning about 8 pressure scale heights (albeit at somewhat higher altitudes than where we are 
making our model comparisons). 

We begin with a 70 = 1.5 model with some scattering present (^ = 0.7, A ~ 0.1) and without including the effect of collision-induced 
absorption (e = 1), using equation l l45t . as shown in the left panel of Figure[8] It is clear that while such a model is roughly consistent with the 
increase in temperature from ~ 0.1 bar to ~ 1 mbar, it fails to produce the increase in temperature with pressure at P > 0.1 bar. By contrast, 
a model with no scattering (^=1,^ = 0) and which includes collision-induced absorption (e = 3500) is roughly consistent with both of 
these trends, but does not quite pass through the data points. As we retain e = 3500 and allow scattering to be present (^ = 0.7, A ~ 0.1), 
the lower atmosphere becomes cooler, thus allowing the model to be consistent with the data points. Too much scattering (^ = 0.5, A ~ 0.2) 
results in a temperature-pressure profile which is cooler than indicated by the data points. 

It is reasonable to ask if our model comparisons are sensitive to the assumed values of the opacities. In the right panel of Figure[8] we 
retain 70 — 1.5 but calculate two more models where we reduce and increase both of the shortwave and longwave opacities by a constant 
factor of 2. It is clear that even with 70 kept constant, the "turn" in the temperature-pressure profile — i.e., first decreasing, then increasing, 
temperature with increasing altitude — is sensitive to the values of the opacities. Lower opacities result in the "turn" residing at higher 
pressures. The temperatures at depth (P > 1 bar) vary sensitively when the opacities are varied due to the changing depths at which the 
stellar photon energy is being deposited (see equation (51). ^^^^^^_^^^^^^^ ^_^ 

For comparison, we include in Figure[8]the temperature-pressure profile for HD 189733b from lMadhusudhan & Seagen ( l2009l) . which 
was obtained via an abundance and temperature retrieval method applied to infrared transit observations from both the HST and the Spitzer 
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Figure 8. Comparison of global-mean temperature-pressure profiles to data points inferred from observations of HD 1 89733b (see text for details and caveats). 
The observed system parameters allow us to set T\^^ = 1697 K, while we set Tint = K and Keg = cm^ g^^ for simplicity. The left panel fixes the 
shortwave opacity and longwave opacity normalization and illustrates the effects due to the variation of the other two parameters. The right panel fixes e and ^ 
and varies the opacities. In the left panel, the horizontal, dash-triple-dot line shows the approximate position of the longwave photosphere, while the hatched, 
yellow region indicates the approximate location of the shortwave photon deposition depth (corresponding to the pair of ^ values adopted). We terminate all 
of the models at P = 0. 1 mbar as our formalism does not treat the physics of the thermosphere. For comparison, we include the temperature-pressure profile 
from Madhusudhan & Seager (2009), labelled "MS09", which was obtained via an abundance and temperature retrieval method. 



Space Telescope. Given the uncertainties associated witli both theory and observation, the profile is broadly consistent with our models for 
P > 0.1 bar but does not include the temperature invers ion from ~ 0.1 bar t o ^ 1 mbar. This is unsurprising because the data point at about 
1 mbar was not included in the original analysis of Mad husudhan & Seageii 1 120091) . 

One may find the status of HD 189733b as the prototypical example of an inversion-less hot Jupiter to be at odds with the increase 
in temperature from ~ 0.1 bar to ~ 1 mbar and our interpretation of it using a 70 > 1 model. This dilemma is resolved if one realizes 
that there are two types of temperature inversions: deep- versus meso-atmospheric inversions. The former refers to temperature inversions 
probed by near infrared water lines, typically at ~ 0.1-1 bar: at these pressure levels, there is general consensus that HD 189733b and HD 
209458b are the pr ototypes o f hot Jupiters w i thout and w ith temperatur e inv ersions in their atmospheres, respectively ( .Burrows et alj|2007 ; 
Charbonneau et al. 2008; For tnev et al.ll2008l : [Madhusud han & Seager 20091) . The latter refers to temp erature inversions from ~ 0.1 bar to 
~ 1 mbar, as probed by the HST transmission spectra of IPont et al.l 1200a) and lSing et al.l d2009ll201 ih . 

Collectively, our models make a crude prediction for the Bond albedo ofHD 189733b: A ~ 0. 1. The Bond albedo is not an observable 
quant ity; rather, it is the geometric albedo which is measured, as has been done for HD 209458b bv lRowe et alj yOOSJ). (See 



Burrows et al. 



20081 for the implications of albedo measurements on modelling.) The conversi on from the geometric to the Bond albedo requires integrating 
over both frequency and phase angle. For example, iBuenzli & Schmid ( l2009i) estimate that the spherical albedq^ As and the geometric 
albedo Ag are related by Ag ~ 0.8.43 in their models of exoplanetary atmospheres with Rayleigh scattering. The reader is again referred to 
§3.4 of ISeageg J2010l) for a detailed discussion of exoplanetary albedos. 



5.2 Semi-analytical models as a guide for three-dimensional simulations 

The ability of our models to predict the temperature at depth (Too) for a given set of parameters provides a useful guide to three - 
dimensional simulations of atmospheric circulation which utilize dual-band, two-stream radiative transfer dHeng. Frierson & Phillippsl201 Ih . 
Such simulations require the specification of the initial temperature and velocity fields. The simplest assumption is to initiate them from a state 
of windless isothermality, which then requires the specification of a constant, initial temperature Tinit- Computational efficiency is optimized 
when the initial temperature field is as close to radiative equilibrium as possible. Physically, we expect the temperature to approach radiative 
equilibrium at depth (P > 10 bar). Therefore, selecting Tinit = Too is a good choice. As the simulation proceeds, the upper atmosphere 
(P ^ 10 b^) adjusts itself to temperatures consistent with dynamical-radiative equilibrium, while the lower atmosphere (P > 10 bar) 
remains in radiative equilibrium. A state of quasi-equilibrium is reached when there is dynamical-radiative equilibrium within the upper 
atmosphere, while the lower atmosphere ceases to transfer significant amounts of energy upwards. Semi-analytical models such as the ones 
presented in this study alleviate the need to perform a tedious parameter search for the appropriate value of Tinit to adopt and therefore 
increase the efficiency of utilizing general circulation models to study the atmospheres of hot Jupiters. 



The Bond albedo is obtained by integrating the spherical albedo over frequency. 



© 2011 RAS, MNRAS 000,[T]tT8 



Clouds and hazes in hot Jovian atmospheres 1 5 



5.3 Future work 



A useful extension of our models will be to consider the dependence of the radius, at a given wavelength, of a hot Jupiter on various 
quantities, 

R = R{T^,Tint,g). (55) 

The temperature at depth Too is a reasonable proxy for the atmospheric effects due to the variation of absorption and scattering. Our formalism 
can then be used to estimate a value for Too, which will provide a lower boundary condition for models of hot Jovian interiors. While the 
elucidation of equ ation (I55t is beyond the scope of the present study, such an exercise will be useful for interpreting the inflated radii of some 



hot Jupiters fe.g., lLaughlin. Crismani & AdamjuOllI) . Specifically, combining our atmospheric models with those of hot Jovian interiors 



will allow a connection to the observed (transit) radius, thereby allowing us to rule out models which over-inflate a given hot Jupiter. 

5.4 Summary 

The salient points of our study may be summarized as follows: 

• We have presented a semi-analytical model for the global-mean temperature-pressure profile of a hot Jovian atmosphere which requires 
the specification of 7 parameters: the intrinsic heat flux (via Tnt), the flux of stellar irradiation (via Trr or Tcq), the surface gravity of the 
exoplanet (g), the strength of shortwave scattering (^ or via A), the shortwave opacity (ks) and the longwave opacity (via kq and e). Allowing 
for an extra source of longwave absortion due to the presence of a uniform cloud/haze layer adds another parameter (kcq), while a purely 
absorbing Gaussian cloud/haze deck adds another 2 parameters to the system (kcq, Ac and Pc; but sets ^ — 1). 

• If a cloud/haze layer contributes an extra source of longwave opacity, then the main effect is to cool and heat the upper and lower 
atmosphere, respectively, analogous to the greenhouse effect on Earth. 

• If a cloud/haze layer causes scattering in the shortwave, then the main effect is to heat and cool the upper and lower atmosphere, 
respectively, thus counteracting its absorbing effect in the longwave — the anti-greenhouse effect. Scattering shifts both the shortwave 
photosphere and the photon deposition depth to higher altitudes. 

• Besides their utility in developing physical intuition, our semi-analytical models provide a guide to three-dimensional simulations of 
atmospheric circulation. In particular, they pin down the temperature at depth for a given set of parameters, which serves as the initial 
condition for the temperature field. 

• We have made a tentative comparison of 4-parameter models to the temperature-pressure data points inferred from the Hubble and 
Spitzer observations of HD 189733b and estimate that its Bond albedo is approximately 0.1. Future observations of HD 189733b will 
corroborate or refute some of our models. 

• The simplicity and versatility of our semi-analytical models allow for an easy comparison to observations. Observers wishing to use our 
models should refer to equation l l45b and set Tnt = K and 7^7^ = as a starting point, thereby distilling the model down to having only 4 
parameters. 
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by a STFC Advanced Fellowship. We thank Isabelle Baraffe, Tristan Guillot, Esther Buenzli, Adam Burrows, Hans Martin Schmid, Sascha 
Quanz, Giovanna Tinetti and Nikku Madhusudhan for useful conversations. 



APPENDIX A: NOTES ON COLLIMATED BEAM AND TWO-STREAM APPROXIMATIONS 

In the coUimated beam approximation, the shortwave intensity may be written as 

Is ifi') = /+ S {fi' ^^i)+I-S {^' + fi) , (A1) 

such that its moments are 

Js = lJ\ Is dp' = i (/+ + /_), 






(A2) 

■^ J— 1 ^ 

1 fl .„ „2 



2, 

where /+ = I+{m) and /_ = I -{in) denote the outgoing and incoming shortwave intensities, respectively. The cosine of the angle 
between the beam and the vertical axis is represented by p, while the Dirac delta function(al), (5(/i' ± /x), describes the sharply peaked 
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Table 1. Table of commonly used symbols 





gcm- 


2 






gcm~ 


2 




bai- 


or dyne 


cm^ 


-2 


bar 


or dyne 


cm" 


-2 


bar 


or dyne 


cm" 


-2 




cm^g- 








cm^g- 








cm^ g- 








cm^g- 








fm2 o- 







Symbol Meaning Units 

cosine of latitude — 

column mass per unit area of atmosphere 

column mass per unit area at bottom of model atmosphere 

vertical pressure 

pressure at bottom of model atmosphere 

pressure level where cloud/haze deck is located 

shortwave scattering parameter (ratio of absorption to total opacity) 

factor for enhancement of longwave optical depth (at P = Pq) due to collision-induced absorption 

shortwave opacity 

longwave opacity 

longwave opacity normalization 

additional longwave opacity due to cloud/haze deck 

longwave opacity normalization of cloud/haze deck cm^ g~ 

ratio of shortwave to longwave opacities — 

ratio of shortwave opacity to longwave opacity normalization — 

ratio of shortwave opacity to longwave opacity normalization due to cloud/haze — 

longwave optical depth — 

shortwave optical depth — 

longwave optical depth normalization in absence of colHsion-induced absorption — 

additional longwave optical depth due to cloud/haze deck — 

thickness parameter of cloud/haze deck — 

blackbody-equivalent temperature associated with internal heat flux K 

irradiation temperature K 

equiUbrium temperature of hot Jupiter K 

global-mean temperature of atmosphere K 

temperature at depth K 

Note: the term "opacity" refers only to absorption opacities. 



angular distributions. A direct consequence of equation | |A2I ( is that Ks/Js = fJ,^ as stated in equation ([9}. Using equations 1 II6I 1 and ( IA2t . 
one may then derive the expression relating the Bond albedo and ^, as stated in equation l IlVI ). by recognizing that 

^=^. (A3) 

For completeness, we examine the longwave intensity in the two-stream approximation by assuming 

Jl (a*') =a + I3fi'. (A4) 

Equation iA4t describes a nearly isotropic radiation field, but in this approximation it is sufficient to describe it by two rays. Taking moments 



A» = 


COS 6 


771 




mo 




P 




Po 




Pc 




e 




e 




KS 




KL 




KQ 




Kc 




KcQ 




■y = 


ks/kl 


70 = 


E Ks/reo 


7c = 


: ks/kcq 


TL = 


- f Kj^ dm 


TS = 


-- Ksm/( 


TO = 


Komo 


T = 


Ki,m 


Tc 




Ac 




Tint 




Ti„ 




-t cq 


= Ti,,/V2 


T 




Tno 





of equation ( IA4t . we get 



fi = ^ = i. £2 = — = ^. (A5) 

Jl 3 ' Jl 3a 



We next derive the ^-values of the two rays in a heuristic fashion; more general derivations exist (e.g. jToon. McKav & Ackermanlll989l) . but 



we are merely seeking to gain physical intuition for the value of the second Eddington coefficient. At the top of the atmosphere, the incoming 
intensity is zero since we are assuming that the star radiates negligibly at long wavelengths. We sample the radiation field (equation IA4I ) 
with two rays, similar to equation l lAU . and denote the outgoing longwave intensity by Jq: 

•Jh - y, -Hl - — , Al - ^~' (.Abj 

whence 

£l=^^\£2 = ^l■ (A7) 

By requiring that £1 — 1/3 be satisfied, we obtain 

m' = i =^ £2 = ^ « 0.58. (A8) 

The term "two-stream" follows from the fact that only two rays at jj. = ±l/\/3 are needed in this approximation. 
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APPENDIX B: NOTES ON THE PHOTON DEPOSITION DEPTH 

Using equations J15t and l l25t . the first moment of the shortwave intensity may be rewritten as 

Hs = ii-oMexp('-^^Y (Bl) 

where Hq = — asBTirr/47r- Taking the latitudinal average, we get 

{Hs) = 2i/o ^3 (\/|ts) . (B2) 

When v^'^'s = 0, we have {Hs{0)} = Ho. We define the photon deposition depth as the layer at which the incident stellar flux diminishes 
to e~^ « 0.368 of its initial value, which is described by 

jSm = 2^^ (^-0 - °-368- (B3) 

If we write the photon deposition optical depth as occurring at \/f rs — vj, then solving equation l lB3t numerically yields vj ^ 0.63. 

APPENDIX C: AN IDENTITY INVOLVING THE DERIVATIVE OF Q 

For an arbitrary function X = X{x, y), we have 

d T"'"^ ,w ^ , T"'"' dX(x,y) , dy2(x)^^, , dy^lx) ^^ , 



We now apply the preceding equation to 

dQ{m,fi,<l)) _ d /""^('"' 



dm 9?n ./„!(,„) 



' (m' , jj.,(j)^ dm'. (C2) 



Since Q — Q{m' , /i, (p) only, we have dQ{m' , /i, cj>)/dm = 0. For a constant m2 — >■ cjo, we have 



dm,o 
lim — ^Q(m2,M,(?!)) = 0. (C3) 

7712^00 dm 



The only non-zero term which remains is 



^ Q{mi,fj.,(t>) = -Q{m,fi,(l>), (C4) 



where we have used mi = m. We thus prove the identity in equation ( 



APPENDIX D: NOTES ON MATHEMATICAL OPERATORS 

As there is some ambiguity in the astrophysical/astronomical community regarding the use of mathematical operators, it is useful to 
clarify the meaning of the operators used in this study. There is generally no confusion over the "=" sign, but when it is used to equate an 
algebraic expression to a numerical answer, it means that no approximations (including rounding-off) are made, i.e., the answer is exact. If 
either an approximation is taken or rounding-off is performed, then the "~" sign is employed. The much-abused "~" sign means "on the 
order of; it does not mean "is proportional to" ("oc"). The "— >" sign means "tends to", i.e., the asymptotic value of a quantity. Finally, the 
"=" sign means "is defined as". 
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